Ising magnets with mobile defects 
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Motivated by recent experiments on cuprates with low-dimensional magnetic interactions, a new 
class of two-dimensional Ising models with short-range interactions and mobile defects is intro- 
duced and studied. The non-magnetic defects form lines, which, as temperature increases, first 
meander and then become unstable. Using Monte Carlo simulations and analytical low- and high- 
temperature considerations, the instability of the defect stripes is monitored for various microscopic 
and thermodynamic quantities in detail for a minimal model, assuming some of the couplings to be 
indefinitely strong. The robustness of the findings against weakening the interactions is discussed 
as well. 



I. INTRODUCTION 

Low-dimensional magnetism in high-T c superconduc- 
tors has attracted much interest, both theoretically and 
experimentallyu. In particular, striped structures in mag- 
nets derived from the LaiCuO^ compound have been 
discussed rather extensively. Motivated by related anal- 
yses and, more specifically, by recent experimentso on 
{Sr,Ca,La)\4 L Cu2AO^i 1 we shall introduce a novel class 
of quite simple two-dimensional Ising models, mimic- 
ing Cu 2+ ions by spin-1/2 Ising variables and holes by 
non-magnetic defects (S= 0). Of course, the aim of the 
present study is not to offer a full, or partial explanation 
of the experimental subtleties. Indeed, beyond the ex- 
perimental motivation, the model is hoped and believed 
to show various intriguing properties being of genuine 
theoretical interest. 

locatWe consider the situation where the spins are 
arranged in chains, with antiferromagnetic interactions, 
J a , between adjacent chains, and a ferromagnetic cou- 
pling, J, between neighbouring spins in the chains, aug- 
mented by an antiferromagnetic coupling, Jo, between 
next-nearest spins in the same chain with a defect in be- 
tween them. The defects are allowed to move through 
the crystal, along the chains. 

The defects tend to form stripes, perpendicular to the 
chains, which, for increasing temperature, first meander 
and then become unstable. To identify the impact of 
the defect mobility on the stripe instability, a 'minimal 
model' is proposed by assuming indefinitely strong inter- 
actions in the chains, J and Jo- This model is studied an- 
alytically, at low and high temperatures, and, for a wide 
range of temperatures, by using standard Monte Carlo 
techniques^. Contact will be made to well known de- 
scriptions of wall instabilities in two dimensions, as have 
been put forward, for instance, in the context-of incom- 
mensurate superstructures in two dimensionsQ and step 
roughening on vicinal surfacesO. Deviations from these 
standard scenarios will be discussed. 

To study the robustness of the properties of the mini- 
mal model and to identify other possibly interesting as- 
pects of this class of Ising models as well, we also consid- 



ered cases with finite couplings in the chains. In addition, 
both for the minimal model and the 'full model', the ef- 
fect of an external magnetic field has been investigated. 

The paper is organized accordingly. In the next sec- 
tion, we introduce the model and elucidate its experi- 
mental background. Then, we present our results on the 
minimal model, followed by a discussion on properties of 
the full model. A short summary concludes the article. 



II. MODELS AND METHODS 

We consider Ising models on a square lattice, setting 
the lattice constant equal to one. Each lattice site (i, j) is 
occupied either by a spin, Sij = ±1, or by a defect corre- 
sponding to spin zero, Sij = 0. The defects are assumed 
to be mobile along one of the axes of the lattice, which 
will be called the chain direction in the following. Neigh- 
bouring sites of the same chain, and (i ± are 
not allowed to be occupied both by defects (short range 
repulsion between defects). The mobility of a defect may 
be influenced by a pinning potential, which, however, will 
be disregarded in most of the following analysis. Even 
for vanishing pinning, the defect cannot diffuse freely, 
in general, because an elementary move, characterized 
by exchanging a defect and a, possibly flipped spin at 
neighbouring sites, is affected by the magnetic interac- 
tions along and perpendicular to the chain direction. We 
assume a ferromagnetic coupling, J > 0, between neigh- 
bouring spins, Si j and Si±ij, along the chain augmented 
by an antiferromagnetic interaction, Jo < 0, between 
those next-nearest spins of the same chain, which are 
separated by a defect. Spins in adjacent chains, Sij and 
Sij±i, are coupled antiferromagnetically, J a < 0. Ac- 
cordingly, the Hamiltonian of the model may be written 
as 



H — — JSijSi±i j + JoSijSi±2j(l — Sf±ij) 

+JgSi t jSi t j±i) — Hy^Sjj (i) 
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where we included a field term. Note that the defects, 
Sij = 0, are separated, along the chains, by at least one 
spin. We shall assume that the number of defects is the 
same in each chain, determined by the defect concentra- 
tion 0, denoting the total number of defects divided by 
the total number of sites. 

The model describes, among others, the thermal distri- 
bution of defects, leading to rather intriguing properties, 
as will be shown below. In particular, at low temper- 
atures, the defects tend to form stripes perpendicular 
to the chains which become unstable at higher tempera- 
tures. The model is believed to be of genuine theoretical 
interest. 

The theoretical model may be motivated by rather re- 
cent experimental findings for the so called telephone 
number compound (Sr, C a, La) 14CU24O 41 which con- 
tains two magnetically one-dimensional elements. One 
subsystem is a sheet like arrangement of Cu^Oz two-leg 
ladders, which is irrelevant in the context of the present 
paper. The second subsystem is an array of CuO-2 chains 
formed by edge sharing CuO^ plaquettes. For this bond 
geometry with a Cu — O — Cu bond angle close to 90 de- 
gree the Goodenough-Kanamori- Anderson rules predict 
a ferromagnetic exchange between nearest neighbour Cu 
ions with spin S = |. This is confirmed by neutron 
diffraction studies of the magnetic structure in the or- 
dered state which, in addition, show an antiferromaa- 
netic coupling of the spins perpendicular to the chainsQ. 
To our knowledge the absolute value of the ferromagnetic 
coupling constant J in the undoped chains has not been 
determined yet from inelastic neutron data. A mean field 
treatment of the magnetic susceptibility suggests a cou- 
pling constant of several meVBEl. 

A long range magnetic order of the Cu spins is only 
observed for certain compositions of the telephone num- 
ber compound. For many compositions the chains con- 
tain a large number of hole- like ckawre carriers. These 
holes imply non-magnetic Cu sitesLnH which inhibit the 
formation of a long range ordered magnetic state. For 
example, in the stoichiometric compound SV14CU24O41 
about 60. percent of the Cu sites in the chains are non- 
magnetiai3'U'E3 aftd-the remaining spins form nearly in- 
dependent dimersoty. The analysis of this dimer state 
shows a rather large antiferromagnetic coupling |Jo| ~ 
11 meV between-. the two Cu spins adjacent to a non- 
magnetic Cu siteoOE^I. This coupling is about one order 
of magnitude larger than the antifeppmagnetic coupling 
between Cu ions in adjacent chainsO. 

The experimental data mentioned so far yield at least 
three relevant magnetic coupling constants in the tele- 
phone number compounds: a ferromagnetic coupling J 
of several meV between nearest neighbour Cu ions in the 
chain, the antifetcomagnetic | Jo| — 11 meV for Cu spins 
adjacent to holescia and an antiferromagnetic coupling J a 
between Cu spins in adjacent chains which is of the order 
of 1 meVlLil. For undoped chains the magnetic proper- 
ties depend mainly on J and J a , whereas the behaviour 
for large hole content is determined by a single coupling 



constant, the antiferromagnetic exchange Jo- For small 
hole concentrations all three magnetic interactions and 
their interplay should be relevant. Experimentally such 
a situation is realized in La^Ca^Cu^iOn where the hole 
content in the chains amounts to about 10 percentile. 
Studies of this compound reveal a very unusual suppres- 
sion of the magnetic order in external fields which can not 
be explained in terms of conventional spin modelsa. It is 
tempting to attribute this strange behaviour to a mag- 
netic field induced movement of the charge carriers which 
frustrates the antiferromagnetic interchain coupling. 

As will be shown below our numerical results for the 
model Eq. (1) indeed reveal a movement of holes in ex- 
ternal fields. We mention that the treatment within an 
Ising model is also related to experimental findings for 
the telephone number compound. Different experimen- 
tal data fopi-Lightly doped chains show a strong Ising-like 
anisotropyntfl which was predicted by two independent 
theoretical treatments fox. spin chains formed by edge 
sharing CuO^ plaquettesrJ. 

To study the above Ising model with mobile defects, 
Eq. (1), we applied analytical low and high temperature 
considerations, and Monte Carlo techniques monitoring 
various thermodynamic and microscopic properties. In 
the simulations, we took into account flips of single spins 
as well as hops of a defect to a neighbouring site in the 
chain leaving a spin at the former defect site. Of course, 
simulations are performed on finite lattices with L x M 
sites (L refers to the chain direction). Usually, we em- 
ployed full periodic boundary conditions. In a few se- 
lected cases free boundary conditions perpendicular to 
the chain direction were applied, corroborating the re- 
sults for periodic boundary conditions to be presented in 
the following. To investigate finite size effects, the linear 
dimensions, L and M, were varied from 20 to 160. Typi- 
cally, runs of at least 10 6 Monte Carlo steps per spin were 
performed, averaging then over a few of such realizations 
to estimate error bars. The concentration of defects, 0, 
ranged from zero to 15 percent. In most cases, we set 
6 = 0.1. 

Physical quantities of interest include the specific 
heat, C, determined from energy fluctuations and the 
temperature dependence of the energy, the magnetiza- 
tion per site, m, and the correlation functions paral- 
lel, Gi(r) = (J2(Si.jSi+r,j))/ LM , and perpendicular, 
ij 

Ga(r) = (%2(SijSi t j+ r ))/LM, to the chain direction. We 

ij 

also calculated other microscopic quantities describing 
the stability of the defect stripes and the ordering of the 
spins and defects in the chains. In particular, we com- 
puted the average minimal distance, d m , between each 
defect in chain j, at position and those in the next 

chain, at + 1) (i.e. ^(min|i — dividing this 
sum over the defects by their number), the cluster distri- 
bution, nd{l), denoting the probability of a cluster with 
I spins of equal sign in a chain (in analogy to .the dis- 
tribution of cluster lengths in percolation theoryll3) , and 
the normalized number of sign changes of neighbouring 
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FIG. 1: Typical Monte Carlo equilibrium configurations of 
the minimal model, O = O.land H — 0, of size L = M — 40 
at temperatures fcsT/|J a | = 0.6 (a), 2.6 (b), and 4.0(c). 

spins, n c , in the chains. Finally, it turned out to be quite 
useful to visualize the microscopic spin and defect con- 
figurations as encountered during the simulation. 

In case of the magnets mentioned above, the absolute 
values of both J and J are large compared to the cou- 
pling between chains, J a . To describe then the behaviour 
of the Hamiltonian, Eq. (1), at low temperatures one 
may consider a simplified model in which the spins form 
intact clusters in the chains between two consecutive de- 
fects changing the sign at the defect, i.e. J and | Jo I are 
assumed to be indefinitely strong. Quantities can now be 
expressed in terms of fcsT/| J a \. The thermal excitations 
are due to motion of the defects in the chains. Again, 
defects are separated by at least one spin. The analy- 
sis of this 'minimal model' will be presented in the next 
section. 



III. PROPERTIES OF THE MINIMAL MODEL 
A. Monte Carlo results 

In the ground state of the minimal model with van- 
ishing external field, H = 0, the defects form straight 
lines perpendicular to the chains and separating antifer- 
romagnetic domains of spins. The ground state is highly 
degenerate, with the degeneracy depending on the con- 
centration 8 of defects. Each arrangement of straight 
defect lines, with a separation distance between the lines 
of at least two lattice spacings, has the same, lowest pos- 
sible energy, resulting in a fast decay of the correlations 
Gi parallel to the chains, while the spins are perfectly 
correlated perpendicular to the chains. The degeneracy 
may be (partly) lifted, for instance, by introducing a pin- 
ning potential or by applying an external field as will be 
discussed briefly below. 

Increasing the temperature, T > 0, the defects are al- 
lowed to move so that the stripes start to meander and 
finally break up, as exemplified in typical Monte Carlo 



FIG. 2: Cluster distribution n d (l) at ksT/\J a \— (circles), 
0.5 squares, 1.0 (diamonds), 1.5 (triangles up), 2.5 (triangles 
left), and 4.0 (triangles down), for the minimal model, G = 
0.1, and H = 0, of size L = M — 40. Results have been 
obtained by exact enumerations at zero temperature, and by 
simulations otherwise. 



configurations depicted in Fig. 1 . It seems plausible that 
the destruction of the defect stripes by thermal fluctu- 
ations is accompanied by singular behaviour of thermo- 
dynamic quantities, like the specific heat, and various 
correlations functions. This suggestion is, indeed, sup- 
ported by the numerical evidence discussed below. The 
effect of both phenomena, meandering and breaking up 
of the stripes, on various physical quantities are shown in 
Figs. 2 to 6, for the case = 0.1. Note that in most of 
the figures we did not include error bars being, typically, 
not larger than the size of the symbols. 

At low temperatures, deviations from the straight 
stripesnipay be characterised by kinks and kink-antikink 
pairsulatSP . Actually, the minimal model resembles 
closely a terrace-step-kink (TSK) model describing step 
fluctuations on vicinal surfaces. The energies of the el- 
ementary excitations may be readily calculated. For in- 
stance, a kink with depth of one lattice spacing costs 
— J a , a kink-antikink pair, created by moving a single 
defect by one site away from the perfect stripe, costs 
— 2J a , for a further diffusion of that defect by another 
lattice spacing away from the stripe an additional — 4J a 
is needed, etc. The fact that consecutive defects in a 
chain cannot be closer than the minimum distance of two 
lattice spacings leads-lo the well known phenomenon of 
'entropic repulsion 'EZlEa between meandering neighbour 
stripes. Due to the entropic repulsion, the meandering 
stripes tend to approach their average distance as given 
by the concentration of defects, 9 (here, at <d = 0.1, 
the average distance is ten lattice spacings). This fea- 
ture is seen, e.g., in the thermal behaviour of the clus- 
ter distribution, nd(l), where the distance between two 
neighbouring defects in a chain is, in the minimal model, 
equal to (I + 1). As shown in Fig. 2 for the moderate 
size L = M — 40, at zero temperature, n<i(l), calculated 
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FIG. 3: Parallel correlation function G\(r) at ksT /\J a \= 
(circles), 0.5 (squares), 1.0 (diamonds), 1.5 (triangles up), 
2.0 (triangles left), and 2.5 (triangles down), for the minimal 
model of size L = M = 60, as obtained from exact enumera- 
tion, T = 0, and simulations otherwise. 

numerically by taking into account all degenerate ground 
states, decreases monotonically with the separation dis- 
tance When turning on the temperature, the maximum 
of rid(l) moves towards equidistant spacing of defects, 
here 1 = 9, and the shape of n<i{l) may be approximated 
by a Gaussian or Wigner function, as has been discussed 
recently in the pantext of terrace width distributions for 
vicinal surfacealj. Increasing the temperature further- 
more, T — > oo, the defects take random positions in the 
chains, and na{l) acquires, of course, again the same Pois- 
sonian form as at T = for random distribution of defect 
lines. The destabilization of the stripes is indicated, for 
example, by a rather rapid decrease, with temperature, 
of the probability to find defects at their average spacing, 
i. e. rid(l = 9) for 6 = 0.1, see Fig. 5. 

As depicted in Fig. 3, meandering and breaking up 
of the stripes may be also observed in the correlation 
function parallel to the chains, G±(r). Again, the be- 
haviour at zero temperature has been determined nu- 
merically without difficulty, for fairly small system sizes, 
by averaging over all ground states. The correlations 
are seen to decay rapidly. The oscillations in Gi(r), al- 
ready hardly visible for L = 60, as shown in Fig. 3, be- 
come less and less pronounced when enlarging L at fixed 
8(= 0.1). The asymptotics of G\{r), for large L, may 
be determined analytically, as discussed in the following 
subsection. Raising now the temperature, T > 0, the 
correlations first become stronger, reflecting the ordering 
tendency which favour equidistant stripes due to the en- 
tropic repulsion, and then decrease quite drastically due 
to the thermal destabilization of the stripes. In fact, the 
perpendicular correlations, G%, fall off rather rapidly in 
the same range of temperatures. Finally, when T — ► oo, 
one encounters again the behaviour at zero temperature, 
with the defects at random positions. 

The destruction of the defect stripes is detected di- 



FIG. 4: Distance between defects in adjacent rows d m (T), 
simulating the minimal model of size L = M — 20 (circles), 
40 (squares), 60 (diamonds) and 80 (triangles up). 
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FIG. 5: Temperature and size dependence of probability for 
next-nearest neighbour pairs of defects, nd(l = 1) (open sym- 
bols) and pairs at average distance rid(l = 9) (full symbols), 
simulating systems of size L = M= 20 (circles), 40 (squares), 
and 60 (diamonds). 



rectly in the average minimal distance between defects 
in adjacent chains, d m . Obviously, d m is equal to zero 
at T = 0, and d m w 2 exp(— | J a \/k B T) at low tempera- 
tures, ksT <C \J a \. In Fig. 4, data for various system 
sizes, L = M ranging from 20 to 80, are displayed. While 
at low temperatures, d m (T) does not depend, in fact, sig- 
nificantly on the system size, it starts to rise rapidly at 
some characteristic temperature, with the height of the 
maximum in the temperature derivative of d m increas- 
ing strongly with larger system size. The location of the 
maximum, at T™^, signalling the breaking up of the 
stripes, moves to lower temperatures as L gets larger. 
The quantitative behaviour is quite similar to the one of 
the specific heat, to be discussed below. 

One possible reason for the destabilization of the 
stripes are effectively attractive interactions between 
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FIG. 6: Specific heat, C, for systems of size L = M — 20 
(circles), 40 (squares), and 80 (diamonds). 



neighbouring defects or lines, mediated by the spins. In- 
deed, such an interaction may occur, for instance, for 
strongly fluctuating stripes so that three consecutive de- 
fects in one chain, j, are in the cage formed by pairs of 
defects in the adjacent chains, j ± 1. Consequently, two 
of the three defects tend to form a pair of next-nearest 
neighbouring defects, as may be checked easily. In any 
event, the probability of such pairs of defects is obviously 
given by nd(l = 1). Its temperature and size dependence 
is depicted in Fig. 5 (together with that of nd{l = 9), 
as mentioned above), showing a drastic increase close to 
the characteristic temperature of the breaking up of the 
stripes, T™ ax . Note that this type of stripe instability is 
not included in the standar d . de scriptions of wall insta- 
bilities in two dimcnsionsoDiIaeJ, where either the num- 
ber of walls is not fixed, giving rise to incommensurate 
structures, or dislocations play a crucial role, in the con- 
text of melting of crystals. Also the bunching of steps, 
in TSK models with attractive step-step interactionsEll 
or instabilities in polymer filaments due to attractive 
couplingsE3'E3 are quite different from the loss of stripe 
coherency we observe here. Of course, the breaking up of 
the stripes has to be distinguished from their meandering 
which may result in their roughness, driven by capillary 
wave excitationsO. 

Meandering and destabilization of the stripes also show 
up in the specific heat, C, see Fig. 6 for systems with 
L = M sites, L ranging from 20 to 80. For each size, C 
exhibits two maxima. The maximum at the lower tem- 
perature is almost independent of the system size, and it 
is related to the kink excitations of the stripes. The upper 
maximum, occurring at T™ ax {L), signals the instability 
of the defect stripes. Its height increases with increasing 
system size, indicating possibly a phase transition in the 
thermodynamic limit, L — ► oo. To estimate the tran- 
sition temperature, we plotted T™ ax versus 1/L, with L 
going up to 160, see Fig. 7. From a linear extrapola- 
tion we obtain approximately fcBT™ ra (L = oo)/|J n | = 



FIG. 7: Size dependence of the location of the maximum in 
the specific heat, T™ ax {L), as obtained from simulations of 
the minimal model, 6 = 0.1 and H = 0, for L = M = 20, 40, 
60, 80, 100, and 160. 

1.05 ± 0.05. Note that finite size analyses, usually for 
sizes up to 80, for other quantities, like G\, G2, and 
rid(l = 9), lead to similar estimates for the possible tran- 
sition temperature. However, the close agreement may be 
fortuituous, depending on the type of the transition. For 
instance, for a Kosterlitz-Thouless transition, the peak 
in the specific heat does not occur exactly at the tran- 
sition temperature, as L — ► 00. A detailed analysis of 
this subtle feature is, however, beyond the scope of the 
present study. 

From simulations of the minimal model with 9 = 0.05 
and = 0.15, L = M = 40, we infer that the character- 
istic temperature, at which the stripes become unstable, 
gets smaller when the concentration of defects, O, is in- 
creased. This observation may be explained by the fact 
that the effectively attractive interactions, caused, for in- 
stance, by the cage effect described above, may set in ob- 
viously at lower temperature when the average distance 
between the stripes decreases. 

Applying an external field, H > 0, the ground states 
change when H exceeds specific critical values. For 
\2J a \ > H > |J Q |, the stripes are no longer straight, 
but they form a zig-zag structure. In that structure, 
supposing the field favours the '+' spins, each '+' clus- 
ter comprises two more spins than the '— ' clusters di- 
rectly below and above that cluster in the two adjacent 
chains. Defects bounding these '— ' clusters are located 
exactly below and above the first and last spins of the 
'+' cluster. Spins and defects in each second chain are 
arranged identically. Obviously, the zig-zag structures 
carry a non-vanishing net magnetization. The degen- 
eracy of the ground state is still high, albeit somewhat 
smaller than in the case of straight stripes at H < |J |, 
because the minimum length of '+' clusters is now three, 
instead of one. For larger fields, H > 2\J a \, the '— ' clus- 
ters shrink drastically: '— ' spins occur only in the pairs of 
next-nearest neighbouring defects; all other spins point 
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in the direction of the external field. 

Monitoring various quantities at fixed defect concen- 
tration, 9 = 0.1, the stripes are observed to become 
unstable at lower temperatures for stronger fields, < 
H < 2\J a \. Varying the field at fixed small tempera- 
ture, for the same defect concentration, we found, that 
the destruction of the stripes seems to be accompa- 
nied by a fairly rapid increase in the magnetization, 
m, leading to an anomaly in the field derivative of 
the magnetization, similar to experimental findings on 
(Sr, Ca, La)nCu2AOi^u. The change from straight to 
zig-zag stripes, on the other hand, leads to a jump in 
m{H) at H — \ J a \ and T = 0. The corresponding maxi- 
mum in the field derivative of the magnetization is, how- 
ever, extremely weak already at very low temperatures. 
More detailed investigations are needed to clarify the ex- 
perimental relevance of these observations. They are be- 
yond our present scope. 

The high degeneracy of the ground states may be lifted 
by introducing a pinning potential. For example, in re- 
lated simulations for H = 0, we found that meander- 
ing and breaking up of the stripes seem to be quali- 
tatively not affected by a weak one-dimensional, along 
the chain direction, harmonic regular pinning potential. 
More realistically, one may introduce a random pinning 
potential. If it is sufficiently strong, it is expected to 
destroy the defect lines even at zero temperature. In a 
weak random potential, the coherency of the defect lines 
gets lost or the lines-jcollide only on a large spatial scale, 
the 'Larkin lengtho. At some finite temperatures, the 
collision length between neighbouring lines due to ther- 
mal fluctuations becomes smaller than the Larkin length. 
Starting from this temperature, the random pinning po- 
tential can be neglected in thermodynamics, though, it 
can be dominant for dynamic phenomena. Thus, our 
model is thermodynamically robust with respect to weak 
random potential except of a very low temperature re- 
gion, in which probably a glassy state occurs. 

B. Analytical results 

Here we derive the asymptotics of the spin correlation 
functions in the chains, G\{r) = (Si_jSi+ r j) , for suffi- 
ciently large distances r, first at zero temperature, also 
valid at infinite temperature, and then at non-vanishing, 
but small temperatures. We start with the obvious state- 
ment that 

SijS^j = (-!)"« (M+r) (2) 

where no(i, i + r) is the number of defects or 'zeros', in 
the interval (i, i + r). Eq. (2) can be rewritten as 

Si,jS i+rJ = (e^MM+r) + e -i^n (i,i+r)y 2 ^ 

Assuming that (rio(i,i + r)) = fio(r) is sufficiently 
large, we apply the Gaussian statistics to the deviation 



<5no(r) = no(i,i + r) — no(r). Thus, the asymptotic ex- 
pression of the correlation function reads 

Gi(r) » cos7m (r)exp (- y ((5n (r)) 2 )^ (4) 

The average uq (r) is related to the concentration of ze- 
ros Co = 6 by fio(r) = cqt. The calculation of ((6na(r)) ) 
is more tricky. First we calculate, at zero temperature, 
the probability p(no,r) of a fixed value of uq at fixed 
r. We find, disregarding here the value of the minimal 
distance between defects, 

p(no,r) = exp(-r(cln— + (1 - c) In - — — ))/ 
c 1 — c 

y/2nrc{l - c) (5) 

with c = no/r. For large numbers, the probability (|^) 
has the Gaussian form near the maximum leading to the 
final result: 

((n - n {r) 2 ) = rc (l - c ) (6) 

Plugging this expression into Eq. (4), we obtain the 
asymptotics of G\{r) = (Si t jSi+ r ,j) a t zero temperature 



Gi(r,T= 0) « cos(7rc r) exp ^yc (l - c )rj (7) 

The asymptotics is valid for r 3> r c , with the correla- 
tion length r c = 2/(ir 2 co(l — Co). The condition of having 
a minimal distance of q + 1 spacings between consecutive 
defects can be taken into account by replacing C( ) by 
C(o)/(l — C(o)<?). For our model q = 1. 

Note that the above considerations hold also in the 
high-temperature limit, T — ► oo, as mentioned before. 

At finite, but small temperatures, fc^T <C \J a \, the 
long-distance asymptotics of the correlation function 
G\ (r) changes dramatically due to the meandering of the 
defect stripes. This process may be described by the free- 
fcrmion approximatiorcl. In this approach the lines are 
represented as trajectories of free fermions. Their en- 
tropic repulsion is treated as statistical repulsion of the 
fermions. The meandering of the lines means that, going 
from one moment of discrete time to the next one, each 
fermion can move to the left or right by one site with the 
amplitude (or probability) z = exp(— 2\ J a \/ksT). These 
processes are described by the Hamiltonian 

H = -z^2(af +1 ai + a+a l+1 ) (8) 

i 

where the fermion operators <Zj and af obey cyclic 
bondary conditions a\^ L = a^ + \ The Hamiltonian is di- 
agonalized by Fourier transformation 
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H = — 2 z (cos p)a+a p ; 
v 

p = 2irm/L(m = 1,2, ...L) 



where 



a k 



E 



e lpk a, 



(9) 



(10) 



The energy band in p-space extends from — tt to tt, 
but it is filled only partly, from — pp to pp, where pp = 
ttco (to take into account the specific minimal distance 
between defects, Co has to be substituted here and in the 
following as before). It means that in the ground state 
(a+ a p ) = 1 for |p| = pp, and (ai a p )= for pp < \p\ < 
tt. In this approach the number of zeros between i and 
i + r is given in terms of the fermion operators by 



n (i,i- 



2+r 

k—i 



(ii) 



The average of this quantity is obviously equal to rco , as 
at zero temperature, but its variance, (Sn 2 ) = ((n (r) — 
no(r)) 2 ), is drastically different. Indeed, 



(Sn 2 ) = ^2 (■ ata-k ■■ at' a W '■) 



(12) 



k.k'—i 



where :XY: denotes the normal product of the operators 
X and Y. Applying the Wick theorem, one obtains 



i+r 



(5n 2 ) = { a k a 'k)( a kat') 



(13) 



fe,fe' 



with the same summation limits. The diagonal term of 
this sum, k = k', gives the contribution rco(l — Co), as at 
zero temperature. However, it will be completely com- 
pensated by the non-diagonal terms. Indeed, the simuL 
taneous correlation function for free fermions is knownO 
to be equal to 



(a k a+) = sin(p F (fc - k'))/(-K(k - k')) (14) 



It is easy to check that (aka^,) 



(ala k , 



In the 



limit of small pp, the summation in Eq.(13) can be re- 
placed by an integration, which can be explicitly per- 
formed without difficulty under the additonal condition 
ppr > 1, giving 



(Sn ) nondiag 



-ppr/n + (lnr)/47r 2 



(15) 



The first term compensates the diagonal contribution. 
Thus, at < ksT ^ |J a |, the correlation function 
Gi(r) = (SijSi+ r j) is finally approximated as 



Gi(r,T= 0) w cos(7rc r)/Vr 



(16) 



Of course, the algebraic decay of the correlations is 
in accordance with-previous findings on free fermions in 
two dimensionsoUiiij. The presence of a pinning potential 
is expected to establish long-range order in the correla- 
tions. Formally, there is no continuous change from the 
exponential decay of the correlations, at T — 0, to the al- 
gebraic decay at non-vanishing temperatures. However, 
the final expression, Eq. (16), is valid only for a sys- 
tem whose size, L, perpendicular to the stripes exceeds 
the collision length of the fermions, l co u = 1/(zcq), which 
goes to infinity as T — > 0. At a fixed value of L, there ex- 
ists a crossover temperature T cr « | J a |/ln(Lcg) at which 
the exponential decay of the correlations goes over into 
an algebraic one. Note that the size, L 3> 1/(zcq), is 
rather large at the concentration we mostly considered 
in the simulations, cq = = 0.1, and the crossover effect 
plays no role there. However, such sizes are not big in 
experimental systems. 

At large temperatures, the correlations G\ are believed 
to decay exponentially, see Eq. (7) which also holds at 
infinite temperature. Thence one expects a transition 
from algebraic decay at low temperatures to an exponen- 
tial decay at high temperatures, in accordance with the 
simulational results. 



IV. BEYOND THE MINIMAL MODEL 

In the following, we shall present results on the full 
model, Eq. (1), with finite ferromagnetic couplings, J, 
between neighbouring spins in a chain and antiferromag- 
netic couplings, Jq, between next-nearest spins in a chain 
separated by a defect. In the simulations we choose 
Jo/ J = —6.25 and J a /J ranging from zero to minus one. 
This choice is, again, motivated by the experimental find- 
ings mentioned above. 

At zero temperature and small fields, one obtains 
the same highly degenerate ground states of perfectly 
straight or zig-zag stripes as in the minimal model. Like- 
wise, the behaviour at low temperatures and H < 2\J a \ 
is characterised by the meandering of the stripes as in 
the minimal model, followed by the stripe instability at 
higher temperatures. However, the instability may be 
masked, for instance, in the specific heat at J a /J = —0.3, 
O = 0.1, and H = 0, for systems of sizes up to L= 80, as 
depicted in Fig. 8 and to be discussed in the following. 

In this case, in addition to the weak, almost size- 
independent maximum at low temperatures due to 
stripe meandering, the specific heat displays a rather 
pronounced peak at higher temperatures being much 
stronger than the one in the minimal model. The peak, 
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k B T/J 



FIG. 8: Specific heat, C, of the full model, Eq. (1), at J a /J = 
—0.3, H = 0, and 0= 0.1, as obtained from simulations for 
systems of size L — M — 20 (circles), 40 (squares), and 80 
(diamonds) . 

however, gets smaller when the system size increases. In- 
deed, it is non-critical, stemming from energy fluctua- 
tions by breaking bonds, J, between spins in the chains. 
It persists when setting J a = 0, i. e. in the onc- 
dimcnsional limit exhibiting, of course, no phase transi- 
tion at all. In that limit, the lower maximum in C disap- 
pears, because there are no stripes. Actually, the defects 
and their mobility play an important role in breaking the 
bonds betwen spins in the chains. For instance, when a 
defect moves next to a flipped spin, the spin on the other 
side of the defect will be flipped rather easily, costing 
an energy of 2 J in the ferromagnetic coupling. In con- 
trast, in chain without defects, an energy of 4J is needed 
to create the elementary excitation comprising a pair of 
neighbouring spins with opposite signs. 

We also analysed, for the case depicted in Fig. 8, the 
simulational data for the spin correlations, G\ and Gi, 
the cluster distribution, rid{l), and the average distance 
between defects in adjacent chains, d m . The data provide 
strong evidence that the defect stripes become unstable 
at about the same temperature, measured in fcsT/|J a |, 
as in the minimal model. Thence the impact of the spin 
flips on the location of the stripe instability seems to 
be rather small for this choice of parameters. In princi- 
ple, the spin flips may lead to new mechanisms, different 
from the effectively attractive interactions between the 
defects discussed for the minimal model, to destruct the 
coherency of the stripes. To detect the stripe instability 
in the specific heat at J a /J = —0.3, presumably signifi- 
cantly larger systems have to be simulated. As shown in 
Fig. 8, e.g., for L = M= 80 merely a shoulder in C starts 
to develop, at about the temperature where the minimal 
model shows a peak in C, compare with Fig. 6. Note that 
very long runs, especially for large systems, are needed to 
get sufficiently good statistics for the Monte Carlo data. 

Applying a magnetic field, H > 0, the breaking up of 
the stripes is found to shift to lower temperatures, as in 



0.12 




kj/j 



FIG. 9: Cluster distribution, nd(l = 1), in the full model, 
Eq. (1), at H= 0, 0= 0.1,and J a /J = -0.1 (circles), -0.3 
(squares), —0.5 (diamonds), and —1.0 (triangles), as obtained 
from simulations for systems of size L = M= 40, see Fig. 2. 

the minimal model. 

When choosing a smaller, but non- vanishing ratio of 
— J a /J, one gets closer to the minimal model. We did 
simulations for J a / J = —0.1. In fact, there the instability 
of the defect stripes is also indicated by a maximum in 
the specific heat, as in the minimal model, already for 
rather small systems, e.g., L = 40, followed by the large 
non-critical peak due to the spin flips in the chains. 

On the other hand, when weakening J with respect 
to J a , the stripe instability, as indicated by the rapid 
increase of the minimal distance d m , may occur quite 
close to the pronounced, non-critical maximum in C be- 
ing due to the spin flips in the chains. In particular, for 
J a /J = —0.5 and, especially, —1.0, one then observes 
another clearly visible maximum nearby the non-critical 
peak in C already for small and moderate system sizes, 
e.g. for L= 40, due to the stripe instabilty. The location 
of the instability, measured in units of fceT/J, increases 
with increasing ratio —J a /J. Interestingly enough, near 
the instability, the probability of finding pairs of next- 
nearest neighbouring defects, rid(l = 1), now does not 
show any longer the overshooting phenomenon, com- 
pared to complete disorder at T — ► 00, in contrast to 
the situation in the minimal model and for small ratios 
J a /J, see Fig. 9. The stripes become unstable at tem- 
peratures of the order of the ferromagnetic spin coupling 
J, and then the tendency to form pairs of next-nearest 
neighbouring defects is diminished by thermal disorder- 
ing. 

V. SUMMARY 

In this paper a two-dimensional Ising model with de- 
fects being mobile along the chain direction has been in- 
troduced. Albeit the model has been motivated by recent 
experiments on cuprates with low dimensional magnetic 
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interactions, the model is believed to be of genuine the- 
oretical interest as well. 

In particular, based on analytical, asymptotical con- 
siderations at low and high temperatures as well as on 
Monte Carlo techniques, the model is found to desribe 
formation of defect stripes, their thermal meandering 
and, at higher temperatures, their destabilization. 

The meandering and the instability of the stripes is 
discussed in the framework of a minimal model, assum- 
ing infinitely strong couplings between the spins along 
the chains. The instability is signalled by pronounced 
anomalies in spin correlation functions, in the spin clus- 
ter distribution along the chain, in the specific heat, and 
in the minimum distance between defects in neighbour- 
ing chains. The breaking up of the stripes seems to be 
caused by an effectively attractive interaction between 
the defects mediated by the spins. 

The main features of the stripe instability persist when 
replacing the infinite couplings by, presumbaly, experi- 



mentally more realistic values. However, the anomaly in 
the specific may be masked for rather small systems, and 
thermal disordering and spin flips may also reduce the 
pairing tendency of the defects. 

New experimental data on a stripe instability in 
(Sr, Ca, La)i4Cu2404i, together with a more detailed 
discussion on possible theoretical interpretations, will be 
presented elsewhere. 
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